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Abstract 

We review the renormalization group method ap- 
plied to non-equilibrium dynamics by tracing the 
way how the hydrodynamic equations can be de- 
rived as reduced dynamics of the Boltzmann equa- 
tion as a typical example. 



1 Introduction 

The theory of the nonequilibrium dynamics 
may be regarded as a collection of the the- 
ory how to reduce the dynamics of many- 
body systems to ones with fewer variables pQ. 
In fact, the Bogoliubov-Born-Green-Kirkwood- 
Yvon hierarchy PI can be reduced to the time- 
irreversible Boltzmann equation 0, which is given 
solely in terms of the single-particle distribution 
function for dilute gas systems 4 . The derivation 
of Boltzmann equation by Bogoliubov[lj shows 
that the dilute-gas dynamics as a dynamical sys- 
tem with many-degrees of freedom has an attrac- 
tive manifold JS] spanned by the one-particle dis- 
tribution function, which is also an invariant man- 
ifold 5 . Boltzmann equation in turn can be fur- 
ther reduced to the hydrodynamic equation (Eu- 
ler or Navier-Stokes equation [Hj) by a perturba- 
tion theory like Chapman-Enskog method |7j or 
Bogoliubov's method [H [8|. Langevin equation jH] 
which may be time-irreversible can be reduced to 
the time-irreversible Fokker-Planck equation with 
a longer time scale than the scale in Langevin 
equation [H] . 

Two basic ingredients are commonly seen in the 
reduction of dynamics, which are interrelated with 
but relatively independent of each other: (i) The 
reduced dynamics is characterized with a longer 
time scale than that appearing in the original (mi- 



croscopic) evolution equation, (ii) The reduced 
dynamics is described by a time-irreversible equa- 
tion even when the original microscopic equation 
is time-reversible [TUl 111) . The fundamental prob- 
lem in the theory of deriving kinetic or transport 
equations is to clarify the mechanism of and to 
implement the above two basic ingredients. 

In a few years ago, we showed [T^] that the so 
called the renormalization-group method l3J 1141 
I15j . which is formulated in terms of the classi- 
cal theory of envelopes |16j. gives a unified the- 
ory for the reduction of the dynamics and can 
be used to derive the various transport equa- 
tions in a transparent way; we have also elu- 
cidated that the underlying mathematics of the 
reduction is an explicit construction of the in- 
variant manifold |121 I17j . a notion in dynamical 
system [HJ [Sj. 

In this report, we pick up the problem to derive 
the hydrodynamics equations from the Boltzmann 
equation as a typical problem of the reduction ap- 
pearing in the field of non-equilibrium dynamics. 
Some remarks will be also given on applications of 
the method to other problems such as the reduc- 
tion of the dynamics of Fokker-Planck equation 
and also on an extension to the relativistic case. 



2 Derivation of hydrodynamic 
equations from Boltzmann 
equation 

In this section, we apply the RG method for- 
mulated in |16l I17j to derive Euler and Navier- 
Stokes equations, successively from the Boltzmann 
equation [T2]. 
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2.1 Basics of Boltzmann equation 

Boltzmann equation I18j is an evolution equa- 
tion of the one-particle distribution function 
f(r,v,t) in the phase space, and reads 



m +v -^ = I[f] 



(i) 



Here the left-hand side describes the change due to 
the canonical equation of motion while the right- 
hand side the change due to collisions; 



I[f] 



dv\ J dv' J dviw(v v\\v' v[) 



x\f(r,v',t)f(r,v[,t) 



f(r,v, t)f(r,V!,t) 



(2) 



which is called the collision integral. The 
transition probability w{v v\\v' v'i) has the fol- 
lowing symmetry due to the time-reversal in- 
variance of the microscopic equation of mo- 
tion; w(vvi\v' v'i) = w(v' Vi\vvi). Furthermore, 
the invariance under the particle-interchange im- 
plies the following equality; w(vvi\v'v'i) = 
w(vi v\v' t v') = w(v' t v'\vi v). 

We say that the function ip(v) is a collision in- 
variant if it satisfies the following equation; 



dvcp(v)I[f] = 0. 



(3) 



Owing to the conservation of the particle number, 
the total momentum and the total kinetic energy 
during the collisions, the quantities given by lin- 
ear combination of the following five quantities are 
collision invariant: l,v and v 2 . 

For a collision invariant <p(v), we can define the 
density and the current j' as follows; 



dvip(v)f(r,v,t), 
dvv<p(v)f(r,v,t), 



(4) 



which satisfy the continuity or balance equation; 



d t n v + V • j. 



0. 



(5) 



Thus we have formally the hydrodynamic equa- 
tions as the balance equations for the conserva- 
tion of the particle number, total momentum and 
kinetic energy. These equations are, however, for- 
mal ones because the distribution function / is not 
yet solved: The solution obtained from Boltzmann 
equation will give the explicit forms of the internal 
energy, the transport coefficients and so on. 



The H function is defined as follows: 
H(r,t) = J dvf(r,v,t)(lnf(r,v,t)-l). (6) 

For equilibrium states, the H function is equal to 
the entropy 5 with the sign changed. Defining the 
corresponding current by 

J H (r,t) = J dvvf{r,v,t)(lnf(r,v,t)-l), (7) 



one has the balance equation; 

dH 
~dt 



+ V-J H = J dvl[f] Inf. 



(8) 



This shows that when In / is a collision invariant, 
the H function is conserved. In fact, one can show 
that this condition is satisfied when f(r,v,t) is a 
local equilibrium distribution function as given by 
(fT9|) below. 

2.2 Application of the RG method to 
derive hydrodynamic equations 

To make it clear that the following discussion fits 
to the general formulation given in [17], we dis- 
cretize the argument v as v — > Vi^: Discrimi- 
nating the arguments (r,t) and Vi in f(r,Vi,t), 
we indicate subscript i for the distribution 

function; f(r,v i: t) = fi(r,t) = (f(r,t))i. Then 
Boltzmann equation now reads 

9fi .. (9) 



where 

I[f]i = ^2 w ( v i v j\ v kVl)(fkfl ~ fifj){r,t).(W) 

j,k,l 

Now let us consider a situation where the fluid 
motion is slow with long wave-lengths so that 

dfi 



v i ~ = 0(e) 1 



(11) 



where e is a small quantity, |e| < 1. To take into 
account the smallness of e in the following calcu- 
lations formally, let us introduce the scaled coor- 
dinate r defined by f = er and d/dr = e ■ d/dr. 
Then © reads 



dfi * r _, dfi 



(12) 



which has a form to which the perturbation theory 
given in is naturally applicable. 

In accordance with the general formulation 
given in |17| . we first expand the solution as fol- 



lows; Uf,t) = f(°\r,t) + ef^(r,t) + 



Let 
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fi(r,t;to)) be a solution around t ~ to given by a Denning the inner product between cp and if) by 
perturbation theory with fi(f, to) being the initial 



value at t = to; 

fi(f,t = t ;t ) = fi(r,t ). 
We expand fi(f,t;to) as 



(13) one can show that A is self- adjoint; 

(<p,Aij>) = {A<p,i/,). 



(22) 



(23) 



fi(r,t;t ) 



fP(f,t;t ) + ef^(f,t,t 



+ 



One can further show that the five invariants 

m,v,YV 2 span the kernel of ^4JH]; KerA = 

{m,v,Y v2 }- The other eigenvalues are found to 

and the respective initial condition is set up as be negative because one can show 

follows, , . ... 

{<P,A<p) < 0. (24) 

~~ ^Oj^o) — fi (^)^o)) (l — 0, 1, 2...). (15) -^Yg wr ite the projection operator to the kernel 

as P and define Q = 1— P. 



The 0-th order equation reads 



(o) 



at 



(16) 



Applying the general formulation in |17| . one 
can readily obtain the first-order solution, 



Now we are interested in the slow motion which 
may be achieved asymptotically as t — > oo. There- 
fore we take the following stationary solution, 



(i) 



-(t - t )Pu 



8f^ 
dr 



df (0) 
A" 1 Q«.-^-.(25) 



The perturbative solution up to the e order is 
found to be 



(o) 



at 



o. 



(17) 



f{f,t,t ) 



f (0 \f,t,t ) + e[-(t-t )Pv <>f 



(0) 



which is a fixed point of the equation satisfying 

(18) 



(i[~f ( % = o, 



+A~ 1 Qv 



dr 



dr 

(26) 



for arbitrary r. Notice that (|18|) shows that the 

distribution function /^°^ is a function of collision RG equation^)], 
invariants. Such a distribution function is a local 
equilibrium distribution function or Maxwellian; 



Notice the appearance of a secular term. 

If one stops the approximation and apply the 



df/dt 



o\t =t 



0, 



n(r,t ) 



in 



3/2 



x exp 



2TTk B T(f,t 

m\vj - u(f,t )| 2 
2nk B T{r,t ) 



one has 
d~f i0) 



+ ePv 



Of 



(o) 



0. 



(19) 



(27) 



(28) 



dt dr 

This is a master equation from which equations 
governing the time evolution of n(r, t), u(r, t) and 

t(o) 



Here, the local density n, local temperature T, 

local flux u are all dependent on the initial time T(r,t) in f x " ; , i.e., fluid dynamical equations: In 
to and the space coordinate f but independent of fact, taking an inner product between m, mv and 
time t. 

The first-order equation reads 



dt >J 



(o) 



dr 



(20) 



mv jl with this equation, one has fluid dynamical 
equations as given below. Notice that / is now ex- 
plicitly solved and the energy density e, the pres- 
sure tensor and the heat flux Qi are given as 
follows; 



Here the linear operator A is defined by 



e(r,t) 



/'[/ (0) ]/ (1) 



™di_ 

h Q h 

(Af (1) ) 



Jdv^\v-u\*fV>(r,v,t) 



■f y 

~f=r t] 



(i) 



(21) Pij(r,t) = J 



-k B T(r,t), (29) 
dvm(vi — Ui){vj — Uj)f(°\r,v,t) 
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nk B T(r, t)6ij 
P(r,t)6 ij} 



Qi(r,t) 



m 
dv — \v 
2 



0. 



The third term represents dissipations. Taking 
(30|nner products between the first equation and 

the particle number, the velocity and the kinetic 
' 'venergy, one will obtain a fluid dynamical equa- 
^3]iion with dissipations included, i.e., Navier-Stokes 

equation 18 . 



We have defined the pressure P using the equation 
of state for the ideal gas in the second line. There 
is no heat flux because the distribution function 
f( ' in the formulae is the one for the local equi- 
librium. Inserting these formulae into ©, we end 
up with a fluid dynamical equation without dissi- 
pation, i.e., the Euler equation; 



dp 

di 

d{pu, 

at 



+ V • pu = 0, 



+ 



+ 



_d_ 



P 



+ e) + 



_d_ 

dxi 



(Py + e + P)ui 



(32) 
0,(33) 

= ()34) 



We notice that these equations have been obtained 
from the RG equation (|28|). It should be em- 
phasized however that the distribution function 
obtained in the present approximation takes the 
form 



f(r,v,t) 



f (0) (r,v,t) 

dfW(r,v,t) 



Or 



(35) 



which incorporates as a perturbation a distortion 
from the local equilibrium distribution and gives 
rise to dissipations. 

One can proceed to the second order approxi- 
mation straightforwardly and obtain fluid dynam- 
ical equation with dissipations as the RG equa- 
tion. The perturbation equation in the second or- 
der reads 



d_ 
dt 



(2) 



df 



(i) 



dr 



(36) 



Here, we must make an important notice: We 
have actually used the linearized Boltzmann 
equation 18 neglecting the second-order term of 

in the collision integral: It is known that 
the neglected term produces the so called Burnett 
terms which are absent in the usual Navier-Stokes 
equation jSJ. 

Applying the RG method, we have 



dt 



(0) 



+ ePv-Vf 



+e 2 Pv ■ VA^Qv • V/ (0) = 0. (37) 



3 Summary and concluding re- 
marks 

In summary, we have shown that the fluid dy- 
namical limit of the Boltzmann equation can be 
obtained neatly in the RG method as formulated 
in |17j . It would be intriguing to see what kind 
of the dissipative fluid dynamical would emerge 
when the RG method is applied to the relativistic 
kinetic equations |19j. 

The reduction of a kinetic equation to a further 
slower dynamics appear quite often, reflecting the 
hierarchy of the space-time of nature. A typical 
problem HO] in this category is to derive Smolu- 
chowski eauation|21| from Kramers equation [25]. 
As is expected, it is possible to develop a system- 
atic theory for the adiabatic elimination of fast 
variables in Fokker-Planck equations, which ap- 
pear in the theory of the critical dynamics 23 . 

We hope that we can report a development on 
these problems in near future. 
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